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0\ ■ ABSTRACT 

General relativistic calculations are made of neutrino-antineutrino 
^ ■ annihilation into electron-positron pairs near the surface of a neutron star. It is 

^ . found that the efficiency of this process is enhanced over the Newtonian values 

up to a factor of more than 4 in the regime applicable to Type II supernovae 
. and by up to a factor of 30 for collapsing neutron stars. 
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Introduction 



O 

^ I It has long been realized that the reaction u + u e'^ + e , near the surface of a hot 

^ I neutron star, is of considerable importance to type II supernova dynamics. At late times, 
in particular, several percent of the neutrino luminosity {L^, ~ 10^^ ergs/sec) from the hot 
^ ; protoneutron star is deposited into the stellar envelope via e~^e~ and neutrino-lepton 

I scattering. This energy deposition, together with neutrino-baryon capture, significantly 
augments the neutrino heating of the envelope in a successful supernova via the so-called 



"delayed shock mechanism" (|Bethe fc Wilson 1985| , |Bethe 1990|) . The late time heating of 



the envelope is most important for the r-process nucleosynthesis. 

Another case of interest stems from neutron star collapse. This could occur either 
from accretion ( [Bhattacharya fc van den Heuvel 1991 ) directly from a companion or from 



plunging into the interior of a companion ( [bryer, Benz fc Herant 19961 ). A third possibility is 



suggested in recent studies of close neutron star binaries near their last stable orbit (|Wilson| 
et al. 19961, [Mathews fc Wilson 1997| ). This work indicates that neutron stars in such 
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systems can compress, heat and emit a neutrino luminosity as high as 10^^ ergs/second. A 
fraction of these neutrinos should be converted to e+e~ pairs via pV e^e~ . In all of these 
scenarios, understanding the efficiency of this reaction is crucial to modeling the energy 
dynamics and observable emission, such as the formation of gamma-ray bursts. 

Previous calculations of reaction efficiencies near neutron stars have been based upon 
Newtonian gravity (Goodman, Dar & Nussinov 1987, Cooperstein, van den Horn & Baron 
1986); i.e. they assume -C 1 where M is the gravitational mass of the neutron star 
and R is the distance scale. In fact can be ~ 0.4 for supernova calculations and ~ 0.7 
for collapsing neutron stars ( [Mathews fc Wilson 1997D and the effects of gravity cannot be 
ignored. 

In this paper we develop a relativistic analytical model to compute the efficiency 
of uV —>■ e~^e~ in a high-gravity environment such as near a neutron star. We find that 
relativistic enhancements to the Newtonian reaction rate can be large; up to factors of 10. 



Neutrino Annihilation 



Newtonian calculations of the — > e~^e energy deposition rate near a hot neutron 
star were done by Cooperstein, van den Horn fc Baron (1986) and [Goodman, Dar ^ 



Nussinov (1987)[ The energy deposition per unit time, per volume (Goodman et al. 1987) 



IS 

(l{r)= 11 U{Pu, r)fTj{pTj, r){cr|v^ - ^^Tjle^ETj} ^" ^ d^p^d^pTj (1) 

where and /u are number densities in phase space, \y is the neutrino velocity, and a is 
the rest frame cross section. The term in curly brackets in equation (|l]) is Lorentz invariant. 
Thus it can be calculated in the center-of-mass frame to be 

{(t|v^ - Wy\e^ey} = ^ {euEtj - p^ ■ PvC^f , (2) 

with Gl = 5.29 x 10"^^ cm^ MeV'^ and 

D = 1±A sin^ ew + 8 sin^ 9w ■ (3) 

Here sin^ 9\y = 0.23 and the + sign is for i/g^e pairs, while — is for z/^I^^ and ^^-17^ pairs. Note 
that we have assumed that the mass of the electrons is negligible since, for the applications 
of interest, the energy of the neutrinos is greater than 10 MeV. Now we replace p^c = Su^u 
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and d Pi, = -^de^dVly where Oj, is the unit direction vector and dVL^, is a sohd angle 

^ " / 1 •^'^•^^(^'^ + ev)elde^eldeTj , (4) 

where the angular integrations are represented by 

;i - ■ VtTjfdVL^dnTj . (5) 



eM . //( 



Thus the energy and angular dependancies have been separated. The energy integrals can 
then be performed in their local frame. Taking fy = p-(e^'^/'^-^ + for fermions 







oo 

4 







we finally get 



4 = ^^ffJ^^\ kT)Mr) ex r(r)^e(r) 



where T(r) is the temperature measured by a local observer. All calculations carried out 
thus far are in the local frame. In order to complete the evaluation of T(r) and 0(r), 
gravitational redshift and path bending must be included. 

It is also of interest to calculate the radial momentum density p imparted to the e^e~ 
plasma from the uV luminosity. This is related to the energy deposition of Equation (^) by 



c c)(r) 



where 



$p(r) = \ I hi-n,- n^jfir ■ (n, + n^))dnjn^ (lo) 



2 

and f is the unit vector normal to the stellar surface 
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3. Bending of Null Geodesies 



From |Misner, 'I'horne fc Wheeler] , the path of a zero mass particle, a null geodesic, in 



Schwarzchild coordinates is given by 

1 drV l-2M/r _ 1 

^d^J ^ 7^ ~ 62 ' ( ) 

where r is the distance from the origin, is the longitude, and b is the impact parameter. 
Here the neutron star mass M is expressed in units of distance; i.e. ^ = 1. All other 
units are in cgs. We wish to express Equation (|TT]) in terms of the the angle 6 between the 
particle trajectory and the tangent vector to a circular orbit. In terms of local radial and 
longitudinal velocities Vr and this becomes 

-nu- / t (12) 



This expression can be substituted into Equation (^) and simplified to give 

r cos 6r 



I _ 2M 



(13) 



All points along a single orbit will share the same impact parameter. Thus a particle 
emanating from the photosphere at R with a trajectory denoted by On will have a trajectory 
9r at radius r given by 




ll- 


2M 


r 




2M 



cosOn . (14) 



It is interesting to note that this equation implies a minimum photosphere radius, i.e. 
Rmin = 3M, below which a massless particle (neutrino) emitted tangent to the stellar 
surface (^/j = 0) is gravitationally bound. The present discussion is restricted to the domain 
R > 3M. 

We may now evaluate Equation (1). Defining fi = sin 6*, then f2 = 
(yU, a/1 — COS0, a/1 — /i^ sin 0) and d^l = cos 6d6d(j) we get 



Q{r)=A'K'^ I f [l-2^^fij, + filfil+hl- nl)] dfi^ dfiTj , (15) 

J X J X 
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where we define 



(16) 



r / 1 



The resuh is 



27r2 

e(r) = — (l-x)^(a;^ + 4a; + 5) . (17) 



In the Newtonian limit (M — » 0) this function has the same form as that in Goodman, Par 



Nussinov (1987]1 and [Cooperstein, van den Horn fc Baron (1987)1 . In the hmit r ^ R this 



becomes 0(r) ^(7)*, so the efficiency of this process falls off strongly with radius as the 
streaming particle trajectories become more parallel. 

For the momentum density deposition, the integration of Equation (|T0|) is similarly 
performed and the result is 

$ (r) = — (1 - x)^(8 + 17a; + 12x^ + 3x^) . (18) 
6 

Thus at the surface of the star (r = R) we have = | and the deposited e+e^ plasma 
moves with a high momentum. At large distances, as r cxo, then — ^ 1, as expected. 



3.1. Gravitational Redshift 

The temperature T{r) in Equation (^) is the neutrino temperature at radius r. This 
must be put in terms of an observable quantity such as the observed luminosity L^. To 
do this we start by expressing the temperature of the free streaming neutrinos at radius 
r in terms of their temperature at the neutrinosphere radius R using the appropriate 
graviational redshift. Temperature, like energy, varies linearly with redshift 



1- 


2M 




R 




2M 




r 



Tir) = ^ ^—4, TiR) . (19) 
The luminosity varies quadratically with redshift 
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L^={l-^)LiR), (20) 
and at the neutrinosphere, the local neutrino luminosity for a single neutrino species is 



L{R) = L, + U = 4rf ^^T(i?)4 . (21) 



Combining these, Equation (R|) becomes 



Ana \ 9/2 (22) 



Equation (^) describes the idealized energy deposition in e~^e~ pairs from the reaction 
u + V —>■ + e~ aX radius r above a neutron star of radius R and with a neutrino luminosity 
Loo- The momentum density deposition is still described by Equation (P). 



4. Results 

In order to quantify the total e~^e~ pair energy deposition, we define Q as the integral 
of q, given by Equation (|3), over proper volume. Q is not an observable quantity and not 
a Lorentz invariant in general relativity, however it is a measure of the total amount of 
energy converted from neutrinos to e+e^ pairs at all radii. With this quantity we wish to 
measure the total amount of local energy deposited via uV e~^e~ and thus available for 
neutrino-electron scattering or other such energy exchange processes. For this quantity we 
get 
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Q 



Q- 



R 



1 - 



2M 



28GItt^C{5) ( k 



3c5/l6 



'1 _ 2MN5 ' 
yRj 



'I — 2M^3/2' 



3/2 



(23) 



where x is defined by Equation (|T6D and y = ^. This becomes 



Q51 = 1.09 X 10-^ ^(1) DL^fi?-'/^ , (24) 

where Q51 and L51 are in units of 10^^ ergs/sec, Rq is the radius in units of 10 km, and we 
have defined 



For the Newtonian case, ^(0) = 1, and we recover the result found in |Cooperstein, van den 



pdorn fc Baron (1986)[ Note that Equation (|2^) is for only a single neutrino flavor. Hence, 



Equation (^ must be evaluated and L51 must be summed over each neutrino flavor in order 
to get the total energy deposition due to annihilation of all neutrino flavors. 

A key purpose of this paper is to quantify the general relativistic augmentation of 
local e~^e~ pair energy deposition compared to the Newtonian calculations. In Figure (|I]) 
the ratio QcR/QNewt = ^{^) is plotted over a range of stellar radii R/M. It is somewhat 
surprising to find that the uV annihilation is enhanced by almost a factor of 2 at the rather 
large radius R = lOM, increasing to a factor of 4 at i? = 5M. As will be discussed, this is 
the relevant range for type II supernova models. In the more extreme relativistic regime, 
down to -R = 3M, pertinent to collapsing neutron stars, this enhancement rises to nearly a 
factor of 30. 
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stellar Radius J3_ 

M 

Fig. 1. — Ratio of general relativistic energy deposition Qgr to total newtonian energy 
deposition QNewt around a neutron star as a function of neutron star neutrinosphere radius 
a) in the range relevant to supernovae b) and the more extreme relativistic range, down to 
R = 3M, relevant to collapsing neutron stars. 

In order to show this enhancement as a function of radius, Figure (0) plots ^ for 
several stellar masses As expected, the enhancement is strongest near the surface 
of the neutron star. Thus the overall efficiency of this process will be dependent upon 
the conditions near the surface of the neutron star such as scale-height of the baryonic 
atmosphere. This point will be important in the discussion of supernovae. 
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Fig. 2. — for various stellar masses Neutrino pair annihilation efficiency drops off 
rapidly with radius. Curves are normalized such that the Newtonian case M = is 1 at 
r = R. 

The enhancements discussed here are the result of two different relativistic effects. The 
ffist is path bending of neutrinos emanating from the neutron star neutrinosphere. The uV 
annihilation cross-section varies like the square of the center-of-mass energy and thus favors 
head-on uV collisions. A simple interpretation of the effect of path bending from Equation 
(0) is that it increases the apparent angular size of the star as seen by a neutrino at 
radius r. This increases the number of near-horizontal neutrinos and thus the probability 
of head-on collisions. 

The second effect is that of gravitational redshift. As the neutrinos rise above the 
surface of the neutron star, they cool. Since uV annihilation depends strongly on energy; 
g oc as in Equation (j^), the annihilation efficiency is strongly reduced by nine powers of 
redshift. However, the energy of neutrinos at the neutrinosphere of the star is not directly 
observable. Only the neutrinos at infinity are measured. Thus q is put in terms of an 
observed luminosity L^o in Equation (|22D . The gravitational field, then, requires a greater 
photosphere luminosity L{R) for a given observed luminosity L^o- For an observed L^o 
there is a net gain in q at the neutrinosphere relative to a Newtonian calculation. It is 
worthwhile to note that this effect is important to all luminosity dependant processes near 
a neutron star, not just uV annihilation. 



- 10 - 



5. Discussion 



The simple analytical arguments presented here are meant to illustrate and quantify 
the importance of general relativity in the physics of emission from compact, massive 
objects. The physics near the surface of a hot (T ~ MeV) neutron star is extremely 
complex, with all possible interactions between baryons, leptons, photons and neutrinos 
playing important roles in energy conversion and transport. Here we have considered only 
one process: + z7 ^ + e~. This is clearly only a small component of a complete model 
of the many processes involved. 

However, the reaction considered here is of particular importance in tapping the 
tremendous neutrino flux emanating from a hot neutron star. Quantifying the dependance 
of this reaction on gravity is crucial to understanding the ultimate energy transport from a 
radiating neutron star. 

The present paper considers two problems in astrophysics for which the efficiency of the 
reaction v + 1/ + e~ near the surface of a hot neutron star is important: supernovae 

and close neutron star binaries. A discussion of each in turn follows. 



5.1. Supernovae 

Figure shows the evolution of the protoneutron star radius during a Type II 
supernova, starting from initial collapse, as calculated by the Wilson-Mayle supernova code 
( |Wilson fc IVlayle 1988| , |lVlayle fc Wilson 19891 ). It is evident that the protoneutron star 



shrinks rapidly for a couple of seconds. This contraction finally halts at a radius of about 
10 km. The density scale-height at the surface of the protoneutron star neutrinosphere is 
indicated in figure (|^) by | ^}°^^ |. As the star contracts and gravity becomes stronger, the 
photosphere becomes more sharply defined. Thus we expect the protoneutron star surface 
to converge to the idealized neutrinosphere model considered here. 



In order to explore this idea we define an angular cutoff factor as in |Cooperstein, vaii 
den Horn fc Baron (1987] 



XNeUr) = z(x-lY{x' + Ax + h) , (26) 



where x is defined by Equation ([16|) with M = 0. XNewt{^) is related to 6(r) of Equation 
(p!7|), but only contains relative neutrino direction information ~ (1 — Q^y ■ fi^)^ and is 
normalized such that x = 1 for an isotropic neutrino distribution. Wilson and Mayle have 
evaluated the cutoff factor, Xw-Mi in their supernova code by solving the Boltzmann 
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equation at various times assuming density, temperature and electron fraction are fixed 
functions of radius. Figure shows a comparison between the neutrino beaming from an 
ideahzed Newtonian photosphere and that calculated by the Wilson-Mayle supernova code. 
At early times, 0.2 seconds after the bounce, the baryon extended atmosphere scatters the 
neutrinos above the photosphere, thus rendering the neutrino angular distribution more 
isotropic than is the case for an idealized sharp neutrinosphere. At early times, t < 0.5 
seconds, this effect of scattering overshadows any effect due to bending and thus simple 
neutrino bending is not a valid model for these times. 
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Fig. 3. — The evolution of a I.4IM0 (= 2.1 km) protoneutron star as calculated by the the 
Wilson-Mayle supernova code is characterized by rapid shrinkage to a radius i? 10 km 
and a strong steepening of density profile at the photosphere as indicated by | | . 

After about one second the protoneutron star radius has contracted to less than 20 
km and its surface has become much more sharply defined, with a baryon atmosphere 
scale-height ~ 100 meters. Thus the calculations exhibited here are valid for R ^ 20 
km when the star surface has converged sufficiently toward a sharp boundary and 
neutrino-baryon scattering is minimized. We conclude that after about one second general 
relativistic effects are important to the dynamics of a supernova. In particular, at late 
times, ~ 10 seconds after the bounce, the matter that will ultimately be involved in the 
r-process is being ejected from the protoneutron star. The heating of the ablated matter 
in relativistic wind calculations ( Pardall fc Fuller 1997|) will be greatly enhanced over 
newtonian wind calulations ( [Woosley et al. 1994] ). 




Fig. 4. — Comparison of angular cutoff factor xw-m obtained from the Wilson-Mayle 
supernova code with that of the Newtonian calculation XNewt from Equation (^Bj). As the star 
shrinks through a) radius i? = 50 km at time t = 0.2 seconds after collapse and h) R = 20 
km at t = 0.5 seconds, the surface becomes more sharply defined and better approximated 
ClS db smooth neutrinosphere. 
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5.2. Close Neutron Star Binaries 

Recent numerical calculations of neutron stars in close binary systems near their last 
stable orbit have indicated that the stars may experience a relativistic compressive force 



and subsequent heating ([Wilson et al. 199q , [Mathews fc Wilson 1997| , [Mathews, Marronetti 



fc Wilson 1998| ). It is argued that 0.5 x 10^^ ergs per star of thermal neutrinos might 



be released within a few seconds. This may be a viable engine for gamma-ray bursts 
( [Wilson, Salmonson, fc Mathews 1997 ; 1998). Accurately understanding the efficiency of 



conversion from neutrinos to other forms of matter (i.e. e'^e~ pairs) is crucial to predicting 
the possible gamma-ray burst characteristics. As shown in Figure (|1|), the relativistic 
enhancements of the reaction uV — > e^e~ for a compact neutron star in a binary can be 
large; up to a factor of 30 for R = 3M. The calculations of [Wilson et al. 199^ indicate 
that R ~ 3M is the radius at which the neutrons stars collapse to black holes. Such 
large enhancements of this reaction rate must be taken into account when modeling the 
emission from a hot, compressed neutron star. Since at late time the neutrino luminosity 
rises to ~ 10^^ ergs/second, we estimate that such enhancements will allow more than 10% 
of the energy emitted in neutrinos to be converted to e+e~ pairs. This estimate includes 
electron scattering of the neutrinos by the annihation-produced pairs, which will more than 
double the energy deposition rate. This energy in e^e~ pairs can then go on to generate a 
gamma-ray burst. 



6. Conclusions 

We have shown that the strong gravity near the surface of a neutron star requires a 
relativistic treatment when modeling stellar emission. In particular, when considering the 
reaction uV e^e~ the effects of path bending and redshift can enhance the efficiency by 
up to a factor of 4 for R = 5M relevant for the proto-neutron star at late-times in a type 
II supernovae. A factor of 30 is possible for R = 3M relevant to collapsing neutron stars. 
These calculations on their own are incomplete when considering the total emission from a 
neutron star; one must perform computations including all relevant processes. However, 
these calculations indicate that relativistic effects will enhance the e^e~ pair energy from 
the neutrino luminosity emanating from the neutron star. This added e~^e~ pair energy 
is important for modeling late-time supernovae and for gamma-ray bursts from collapsing 
neutron stars. Thus general relativistic computations are necessary for the accurate 
modeling of neutrino winds from neutron stars. 

Work performed under the auspices of the U.S. Department of Energy by the Lawrence 
Livermore National Laboratory under contract W-7405-ENG-48. 
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